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Abstract. - A number of strange results have been reported for the twist elasticity of a 
semiflexible filament, actin. In particular dynamic and static methods for determining the 
torsional twist modulus give very different results. I show here that internal losses (friction) in 
a semifiexible filament could be an important cause of dissipation and that the interpretation of 
dynamic experiments should be made in terms of a storage and loss modulus for actin proteins. 



Several recent experiments ^ have tried to characterize the torsional persistence length of 
actin filaments. The most straight forward measure the static twist fluctuations by attaching a 
bead to the end of a single filament and following its torsion phase [Q . These experiments give 
a result which is close to that well established for the bending elastic constant measured 
from transverse fluctuations; the result is to be expected with a material with Poisson ratio 
close to zero @. A second class of experiments [D has always seemed to give strange results. 
They depend on measuring the torsional dynamics, and are not in agreement with the static 
measurements. I have been motivated by this anomaly to reconsider the theory of torsional 
dynamics of semiflexible polymers and the importance of dynamic bend-torsion coupling. A 
simple argument will show that there can be non-hydrodynamic contributions to the dissipation 
and that these could be as important as hydrodynamic drag on the twist dynamics. Since 
dynamic methods really measure the product of a stiffness and friction coefficient the results 
of static and dynamic experiments should not necessarily be the same. For this mechanism to 
be important we are forced to admit that actin is described by a viscoelastic modulus with an 
internal mode near IMHz. Since it is well known that many proteins have slow conformational 
changes, this hypothesis appears reasonable, though clearly needs experimental confirmation 
for this particular protein. 

Actin filaments are helical structures with a diameter, d of 7 nm and lengths which can reach 
tens of microns. They have a bending persistence length which is close to 15 microns ^ . Much 
theoretical work has considered the bending dynamics of these filaments however relatively less 
thought has been given to twist dynamics in these systems; however careful theoretical work 
on twist dynamics in DNA ^ has been used to interpret a series of dynamic torsional 
experiments on actin [Q. The twist dynamics can be observed by a technique of polarized 
fluorescent scattering. One marks the sample with a fluorescent probe and illuminates with a 
laser. The dye absorbs a photon, which is slightly later re-emitted. If the dye rotates between 
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the absorption and rccmission of the photon there is a decorrelation in the signal compared 
with a stationary object. 

The rest of this article is organized as follows. We first summarize the rotational dynamics 
of semiflexible polymers in the limit of linear dynamics. We then note that the experimental 
situation in actin is expected to be out of the range of validity of the simple linear theory, in 
particular the dissipation links the twist and bend degrees of freedom. This coupling can lead 
to extra internal losses which we shall estimate in a rather crude manner. 

In linear theory there are two main contributions to the rotation of the dye fixed on a 
filament. Firstly the internal modes of the filament obey a Rouse like equation 

where (f){s) is the angle turned by a filament with respect to the fluid, K is the twist stiffness 
and p0 a friction coefficient, /(s, t) describes the Brownian noise of the fluid. Secondly the 
bending of the filaments in space leads to a rotation of the tangent of the filament. The 
bending dynamics of a filament are described by a Langevin equation |^ 

With K the bending elastic constant, p a dissipation coefficient and rj^ the transverse amplitude 
of the fluctuations of the filament. The direction of the filament is given by its tangent, 
t = dsr{s). 

The rotation of a monomer on a filament (measured in a dynamic experiment) has two 
contributions 

< t) - Hs, 0)r >^ J l-expMXt,Vp.) ,^ ^ (^/^^^)r/. 
from eq. (^) and 

< idsvAt) - dsv^m' f ^2 l-exp(^2^gVp) ^^ ^ {t/pn')'/^ (4) 

J K,q 

from eq. (^) . If we look at the time scales in the eq. (|l|) and (|^) we see that there are two 
relaxation times rt, the twist relaxation time and Tf, the bend relaxation time which scale 
differently with the excitation wavelength, q. In fact r^i/rf ^ (qd)"^ with d the diameter of 
the filament. Experimentally for a ten micron wavelength mode we expect that Tt should be 
several seconds while t* is rather 100/Lts. Given the great disparity in time scales one expects 
that (|3|) is the dominant term in the twist dynamics in most experimental situations (though in 
gels both should be important since the free rotation of the filaments is impossible) . Dynamic 
depolarization experiments are then sensitive to correlation functions of the form 

/(0=exp(-a($(s,i)-$(s,0))2) (5) 

Where $ is a sum of the two contributions (3) and (4) leading to a stretched exponential decay 
of the signal. For more detail see a is of order unity. 

These two independent contributions are derived, however, via linearisation of the dynamic 
equations. We shall now make the argument that for long filaments the twist and bending 
are strongly coupled and that this leads to important modifications of eq. (|l|) and (||) : When 
filaments are very short the simple twist and bend modes are eigenvalues of both the static 
and the dynamic problem, since the friction matrix is diagonal and does not couple the bend 
and twist degrees of freedom. For longer filaments the situation is more complicated. There 
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is no coupling in the energy but there is a dynamic couphng in the friction. As an illustrative 
example consider a crankshaft shaped filament. If the filament is twisted at its ends in a 
resistive fluid a bending force is exerted on the filament. Thus we are force to reconsider the 
nature of the dynamic eigenmodes for long filaments. 

For a long filament one can envisage two extreme situations. Either (a) that the filament 
keeps its shape in body frame and thus changes its shape rapidly in the laboratory frame i.e. 
that it does spin as a crankshaft. Or (b) The filament keeps a nearly constant shape in the 
laboratory frame and spins rapidly in a narrow tube, implying that the filament undergoes 
rapid bending in the body frame. Experimentally one sees under a microscope no evidence 
of case (a), one also understands that rapid rotation of a constant shaped filament will be 
unfavorable due to the high friction coefficient of the spinning object. We conclude that 
superposed on the slow shape fiuctuations of a actin filament is a very rapid spinning of 
the filament which forces the filament to bend at high frequencies, typically 10 kHz. More 
formally we know that a general non-planar curve is characterized by non-zero torsion and 
writhe. Even if we suppress completely the torsional fluctuations (described by the variable 
(/)) geometry implies that bending fluctuations necessarily lead to rotational movement of the 
fllament and thus additional contributions to the dissipation. 

We now argue that because of the large ratio of the relaxation times of the bending and 
twisting modes we can use simple physical arguments to deduce the appropriate dynamical 
variables in this coupled, non-linear problem. We flrst consider the case of a filament embedded 
in a gel (or almost equivalently in a dense solution of filaments). In this case the long 
wavelength bending degrees of freedom are completely frozen out. However the filament can 
continue to twist freely according to eq. if we interpret (j) as rotation in the tube of the gel. 
If we now relax the absolute constraint on the bending degrees of freedom imposed by a gel we 
realize that the large ratio of to Tt leads to an almost identical situation; the bending modes 
can not relax on a time scale of the twist motion. We can deduce that the two dynamic modes 
of the filament in free solution are (1) A rapid twist mode in a narrow tube which is described 
by eq. (|l|), which however leads to strong bending fiuctuations in the body frame; and (2), a 
slow bend mode which is followed adiabatically by the much faster twist degrees of freedom. 
For frequencies lower than a few kHz this mode should be well described by eq. (^. We now 
pass on to the consequences of these remarks on the calculation of the twist dynamics. 

The energetics of semifiexible polymers are usually modeled using a description in terms of 
a simple elastic theory of a uniform rod or more accurately a helical structure Q . In these 
models the filament is taken to perfectly elastic; that is the storage modulus is independent of 
frequency. In general we expect that at high frequencies deformations are going to be couple to 
internal, dissipative modes and that a more accurate description is in terms of a viscoelastic, 
Voigt, solid described by the constitutive equation 

E{lj) = — ^ (6) 

where Eq is the low frequency elastic modulus (measured in static experiments) and r is a 
characteristic internal time scale. Experiments performed at finite frequency are thus coupled 
to both the real and imaginary part of the response function. We shall now show that a simple 
estimate of the dissipation coming from eq. (^) leads to a modification of the friction coefficient 
in (|l|) if the internal time scale in the constitutive equation is slower than 10~^s. We not 
that the case of a broad spectrum of modes should be even more interesting: In particular in 
the case of a material with many low frequency modes the dynamics should become highly 
anomalous. 

The material equation eq. (0) can be written in terms of storage and loss moduli. At 
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frequencies satisfying lot << 1 eq. has the foUowing interpretation: The maximum stored 
energy of an excitation of amphtude x varies as e = Eqx^ /2. However, there is also a weak 
dissipation of energy equal to 

dEi/dt = euj^T (7) 

We can also write the conventional hydrodynamic dissipation in a similar form, so that if we 
spin a rod of diameter d, length I at angular frequency ui in a fluid of viscosity rj the dissipation, 
is equal to 

dE2/dt = TTuj'^d^lr] (8) 

The fact that both Ei and E2 are quadratic in lo shows that we have a renormalization of 
the rotational frictional coefficient for bent filaments. The most delicate part of our argument 
now follows where we shall try to estimate e for a filament which has been bent by thermal 
fluctuations. 

Very short fllaments are essentially straight, they can spin without hindrance along their 
long axis and there is no coupling between the bend and twist modes. Longer filaments are 
described by transverse fiuctuations which scale as 

il - lyip (9) 

with the amplitude of fluctuations, I the length of filament and Ip the persistence length 
[ p] . When these transverse fluctuations have an amplitude comparable to the diameter of the 
filament the rotational friction of filament will increase rapidly with the filament length. For 
filaments longer than this, as we have argued, it is much easier for the filament to bend while 
spinning and the twist eigenmode in the laboratory frame is in fact a combined twist-bend 
mode. Thus the coupling of the bend to the twist will become important when / ^ {d^lp)^/^. 
If we make a filament of this length oscillate an amplitude d we are in turn stocking an energy 
comparable to ksT in each section of length I. Thus we argue that for e we should count 
UbT / {d^lpY^'^ per unit length of filament. We are now in a position to compare the two 
sources of dissipation: The ratio of the two dissipation rates varies as 

These two sources of dissipation become comparable if we let r be equal to 10~^s. As noted 
above this does not appear too unreasonable for a protein, however independent verification 
of this point would be useful. 

A particularly interesting remark in implies that the twist dynamics are only weakly 
dependent on the viscosity of the solvent when this is changed by adding glucose. This would 
seem to be strong evidence that the dissipation is dominated by non-hydrodynamic sources. 
However, the values of the twist stiffness which are deduced dynamically are smaller than 
those deduced in static experiments. From the scaling in eq. (3) of a characteristic time in the 
product p,f,K one would rather expect the dynamic measurement to produce larger values than 
the static experiments. The numerical fitting procedure used in Q is, however, particularly 
complicated including corrections due to short filaments which make a simple re-fitting of 
the data essentially impossible; a full re-analysis of the data should perhaps be done with a 
renormalized value of the effective rotational friction coefficient. 

We note that for length scales smaller than I we should not be coupling the twist and bend 
modes. Thus one can expect some sign of a dynamic crossover for filaments or mode of scale I, 
or at frequencies comparable to kBTlp/l^i] in the bending dynamics. This is probably at some 
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tens of kHz, frequencies which are accessible through DLS or micr orhcology. g |l|. It 
would be interesting to try and make a more quantitative theory of bend twist coupling by 
using the large ratio of time scales as an expansion parameter in a adiabatic expansion of the 
dynamics in order to produce a quantitative theory of this crossover. 

We note that the scale I by numerical coincidence is comparable to the helical pitch of actin 
filaments. One expects that the two principal bending elastic constants of a filament are not 
quite the same due to the low symmetry If the filament is forced to spin in a constant 
shape this could introduce further dynamic couplings which should be investigated. In the 
literature there are also many reports of the presence of relatively rare "defects" or "kinks" 
in actin filaments so that the equilibrium form of a filament consists of long straight segments 
followed by turns by several degrees. Such kinks, if they are present, must have a major effect 
on twist dynamics and can lead to further novel contributions to the relaxation function. The 
presence of such kinks however would considerably complicate the fitting of the experimental 
data. 

We should contrast our discussion here with a recent article in which the twist dynamics 
are hypothesized as being ultra slow due a local conservation law for the writhe. This ultra 
slow twist relaxation in turn couples to the bending modes due to the geometry of twist 
leading to slow relaxation of bending modes. We see here that the torsional modes are very 
rapid compared to the bending modes. This is due to the fact that the relaxation time of the 
bending modes scales as g^, while the twist modes scale more slowly in q^. It seems unlikely 
that such considerations of ultra slow dynamics apply in experiments on actin solutions. 
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